1 



OS 

os 



Extraction of skewed parton distributions from experiment 

Andreas Freund a 
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In this paper we will discuss algorithms for extracting skewed parton distributions (SPD's) from experiment as 
well as the relevant process and experimental observable suitable for the extraction procedure. 
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1. Introduction 

The basic concept of SPD's [ Q, § |, g § is best 
illustrated with the lowest order graph of deeply 
virtual Compton scattering (DVCS) in which a 
quark of momentum fraction X\ leaves the proton 
and is returned to it with x 2 - The two fractions 
not being equal is due to the fact that an on-shell 
photon is produced which necessitates a change 
in the + momentum in going from the virtual 
space-like photon with + momentum —XbjP+, to 
basically zero + momentum of the real 7. This 
sets X2 — x\ — x and thus the skewedness param- 
eter to x. Thus one has a nonzero momentum 
transfer onto the proton and the parton distri- 
butions (PDF's) which enter the process are non 
longer the regular PDF's since the matrix element 
of the quark (gluon) operator is now taken be- 
tween states of unequal momentum rather than 
equal momentum. 

2. Appropriate Process and experimental 
observable 

The most desirable process for extracting 
SPD's is the one with the least theoretical un- 
certainty, the least singular Q 2 behavior and a 
proven factorization formula. 

The process which fulfills all the above criteria 
is DVCS and the experimental observable which 
allows direct access to the SPD's is the azimuthal 
angle asymmetry A of the combined DVCS and 
Bethe- Heftier (BH) differential cross section. A is 
defined as f 01 : 
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The reason why this asymmetry is not is due 
to the interference term between BH and DVCS 
which is proportional to the real part of the DVCS 
amplitude. The factorized expression for the real 
part of the amplitude is [ El 



Re T(x, Q 2 



dv 

-^Re C % {x/y,Q 2 ) 

i+x y 



My,x,Q 2 ). (2) 

Re Ci is the real part of the hard scattering coef- 
ficient (HSC) and are the SPD's. At HERA 
one is mainly restricted to the small- x region 
where gluons dominate and thus i will be g. Thus 
Eq. ([!]) contains only measurable or directly com- 
putable quantities except Eq. (||) in the interfer- 
ence part. Therefore, one would now be able to 
extract the SPD's from on A, if one could decon- 
volve Eq. (||). However, for SPD's life is not as 
"simple" as in inclusive DIS where the deconvo- 
lution for F 2 is trivial, since the SPD's depend on 
two rather than one variable. Furthermore, the 
HSC depends on the same variables as the SPD's. 
These facts make the deconvolution of Eq. (||) im- 
possible. 

3. Algorithms for extracting SPD's 

Rather than deconvoluting, one can expand the 
PDF's with respect to a complete set of orthogo- 
nal polynomials Pj ap \t). In this particular case 
we need the orthogonality of the polynomials to 
be on the interval -1 < t < 1 with t = 
equivalent to— l + a;<y<las found in Eq. 
One can then write the following expansion: 



r> g (y,x,Q 2 ) = 
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M™(x,Q 2 ) 



(3) 



with u>(i|ap), rij(ap) and up being weight, nor- 
malization and a label determined by the choice 
of the orthogonal polynomial used. M 9:9 (x, Q 2 ) 
is given by: 



M]> 9 (x,Q 2 )=J2E% 9 (x)f™(x,Q 2 ), 

fc=0 



j 

where 



f]> 9 (x,Q 2 ) ^^B^f^&Q 2 ). 



(4) 



(5) 
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Bjfc is an operator transformation matrix which 
fixes the NLO corrections to the eigenfunctions 
of the kernels and is thus just the identity matrix 
in LO. The upper limit in Eq. (^]) is given by the 
constraint ^-functions present in the expansion 
coefficients, which are generally defined by 

e jk r»i> + fc) 



E jk (v\a P \x) 
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The moments of the SPD's evolve according to 

ff B {x,Q 2 )=E j {Q\Ql)ff a {x,Ql) (7) 

where the evolution operator is a matrix of func- 
tions in the singlet case. Finally, the Gegenbauer 
moments of the SPD's at Qq are defined by 
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f 9 (t,x,Q 2 ). (8) 



In LO order and at small x the above formal- 
ism simplifies. Owing to the conformal properties 
of the operators involved in the definition of the 
SPD's one finds the following expansion 



f 9 ( yi ,x,Q 2 ) = 



w(t\5/2) 
2-x^^ NAb/2) 



vv 



E] k (x)cf(t)f k (x,Q 2 ), (9) 

with w{t\v) = 0(1 - t)) v - x l 2 and the Cj"s being 
Gegenbauer polynomials . The multiplicatively 
renormalizable moments evolve as above but with 
the explicit evolution operator: 

Ef{Q 2 1 Q 2 ) = Te^ h ^l WMr)) ) (10) 

where T orders the matrices of the regular LO 
anomalous dimensions (i,k = q,g) along the inte- 
gration path. 

Inserting Eq. (^) in Eq. (|J) one obtains: 

oo oo 

Re T(x, Q 2 ) = E k ^(Q 2 , Q 2 ) 

3=0 k=l 

1 dt w(t\h/2) 



f 9 k _ x {x,Ql)E° k _ x {x) 



_ x 2t + x Nj (5/2) 
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(ii) 



where we chose the factorization/renormalization 
scale to be equal to Q 2 . As one can see the inte- 
gral in the sum is now only over known functions 
and will yield, for fixed x, a function of j as will 
also do the expansion coefficients for fixed x. The 
evolution operator can also be evaluated and will 
yield for fixed Q 2 also just a function of j, which 
leaves the coefficients f k _ x (x, Qq) as the only un- 
knowns. Since the lefthand side will be known 
from experiment for fixed x and Q 2 , we are still 
in the unfortunate situation that a number is de- 
termined by the sum over j of an infinite number 
of coefficients. However measuring the real part 
at an infinite number of Q 2 for fixed x, one would 
have an infinite dimensional column vector on the 
lefthand side and on the right hand side one would 
have a square matrix times another column vec- 
tor of coefficients of which the dimension is deter- 
mined by the number of j. Since all the entries in 
the matrix are real and positive definite provided 
that there are no zero eigenvalues one can find 
the inverse. Thus one can directly compute the 
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moments of our initial parton distributions which 
are needed to reconstruct the skewed gluon dis- 
tribution from Eq. (||). 

The drawback of the above procedure is that 
this process has to be repeated anew for each x 
and that NLO evolution studies [ || indicate that 
one might need as amuch as 50 — 100 polynomi- 
als to achieve enough accuracy at small-x. This, 
of course, would render this procedure useless in 
an experimental situation where a j of 5 — 10 is 
possibly achievable! 

A practical way out of the above mentioned 
predicament is by making a simple minded ansatz 
for the skewed gluon distribution in the different 



regions like Aqz Al (l 



for the DGLAP re- 



gion where z is now just a dummy variable, plug 
this form in Eq. (||) and fit the coefficients to 
the data of the real part of the DVCS amplitude 
for fixed x and Q 2 . One can repeat this proce- 
dure for different values of Q 2 and then interpo- 
late between the different coefficients to obtain 
a functional form of the coefficients in Q 2 , al- 
ternatively, after having extracted the values of 
the coefficients for different values of x at the 
same Q 2 , use an evolution programm with the 
ansatz and the fitted coefficients as input and 
check whether one can reproduce tha data for the 
real part at higher Q 2 , thus checking the viability 
of the model ansatz. 

To obtain an ansatz fullfilling the various con- 
straints for SPD's (see Ji's and Radyushkin's ref- 
erences in [ Q ) , one should start from the double 
distributions (DD) (see Redyushkin's references 
in [ ED].) which yield the skewed gluon distribu- 
tion in the various regions 



1-x 



g{y,x) = 6(y>x) dzG(y - xz, z) + 

Jo 

% <x) j dzG(y-xz,z). (12) 
Jo 

Due to the fact that there are no anti-gluons, the 
above formula is enough to cover the whole re- 
gion of interest — 1 + x < yleql. What remains is 
to choose an appropriate model ansatz for G, for 
example, 

h(z\,z) 



with f{z\) being taken from a diagonal 
parametrization with its coefficients now being 
left as variants in the skewed case and the normal- 
ization condition h(zi) = J Q dzh(z\,z) such 
that, in the diagonal limit, the DD just gives the 
diagonal distribution. The choice for h(z\, z) is a 
matter of taste but should be kept as simple as 
possible. 

4. Conclusions and outlook 

After having showed, that the extraction of 
skewed parton distributions from DVCS experi- 
ments is principally as well as practically possi- 
ble. 
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